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Abstract 

Background: Inorganic phosphate is an essential nutrient required by organisms for growth. During phosphate 
starvation, Saccharomyces cerevisiae activates the phosphate signal transduction (PHO) pathway, leading to 
expression of the secreted acid phosphatase, PH05. The fission yeast, Schizosaccharomyces pombe, regulates 
expression of the ScPH05 honnolog [phol^) via a non-orthologous PHO pathway involving genetically identified 
positive iphoT^) and negative (cskl^) regulators. The genes induced by phosphate limitation and the molecular 
mechanism by which phoT^ and cskl^ function are unknown. Here we use a combination of molecular biology, 
expression microarrays, and chromatin immunoprecipitation coupled with high-throughput sequencing (ChlP-Seq) 
to characterize the role of pho7^ and cskl^ in the PHO response. 

Results: We define the set of genes that comprise the initial response to phosphate starvation in 5. pombe. We 
identify a conserved PHO response that contains the ScPHOS iphol^), ScPH084 (SPBC8E4.01c), and ScGITl 
(SPBC1271.09) orthologs. We identify members of the Pho7 regulon and characterize Pho7 binding in response to 
phosphate-limitation and Cskl activity. We demonstrate that activation of phol^ requires Pho7 binding to a UAS in 
the phol^ promoter and that Cskl repression does not regulate Pho7 enrichment. Further, we find that 
Pho7-dependent activation is not limited to phosphate-starvation, as additional environmental stress response 
pathways require phoT^ for maximal induction. 

Conclusions: We provide a global analysis of the transcriptional response to phosphate limitation in S. pombe. Our 
results elucidate the conserved core regulon induced in response to phosphate starvation in this ascomycete 
distantly related to S. cerevisiae and provide a better understanding of flexibility in environmental stress response 
networks. 
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Background 

Inorganic phosphate (Pi) is an essential nutrient required 
for signal transduction, energy metabolism, and biochem- 
istry in all organisms. Maintaining a constant, stable 
concentration of internal inorganic phosphate is a major 
challenge for biological systems. Because external concen- 
trations of inorganic phosphate fluctuate unpredictably. 
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microorganisms have evolved strategies to sense external 
phosphate levels [1-3], communicate this information to 
the nucleus [4,5], and induce transcription to respond to 
phosphate flux [6-8]. The phosphate signal transduction 
(PHO) pathway in the budding yeast, Saccharomyces cere- 
visiae, is the most thoroughly studied example of phos- 
phate homeostasis in eukaryotes [9-11]. 

The transcription factors Pho4 and Pho2 play a key 
role in the phosphate starvation response in S. cerevisiae. 
When cells are grown in conditions where inorganic 
phosphate is plentiful, Pho4 is multiply phosphorylated 
by the cyclin-dependent kinase-cyclin (CDK-cyclin) 
complex, Pho85-Pho80 [12]. When Pho4 is phosphory- 
lated, it is localized to the cytoplasm [13,14], does not 
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interact with Pho2 [15], and the PHO regulon is not 
expressed. During phosphate starvation, the CDK inhibi- 
tor Pho81 binds to the secondary metabolite myo-D- 
inositol heptakisphosphate (IP7) and inhibits the Pho85- 
Pho80 complex [16,17]. Inhibition of Pho85-Pho80 
allows Pho4 to be dephosphorylated, enter the nucleus 
[18], co-operate with Pho2 [7], and induce a set of genes 
responsible for harvesting inorganic phosphate from the 
environment [6]. Pho4 function can be conveniently 
monitored by measuring the activity of the secreted acid 
phosphatase, Pho5, which is one of the most highly 
induced members of the PHO response [19,20]. The 
genes that comprise the PHO regulon have been well 
characterized and the precise sites within the genome 
where Pho4 binds during phosphate starvation are 
known [7,21]. Pho4 regulation occurs in response to 
changes in external phosphate levels and Pho4 activity is 
not thought to be regulated by other stress responses. 

In this study we ask the following: is the PHO tran- 
scriptional response observed in S. cerevisiae conserved 
in the distantly related ascomycete, Schizosaccharomyces 
pombel S, pombe presents an interesting opportunity for 
addressing this question because: (1) 5'. pombe did not 
experience a recent whole-genome duplication event - 
thought to contribute to specialization [22] - possibly 
preventing the PHO response from developing a dedi- 
cated regulatory network; (2) the orthologs for the PHO 
pathway either do not exist {PH08h PH02, PH04) or 
are not involved in the PHO response {PHO80, PH08S) 
in S. pombe [23]; and (3) recent work utilizing a deletion 
collection in S. pombe has outlined a basic regulatory 
structure for the Pi-inducible, secreted acid phosphatase 
phol^ (the ortholog to PHOS) creating an opportunity 
for comparison with the S, cerevisiae PHO response 
[24]. During phosphate starvation S, pombe Pho7 - a 
putative transcription factor containing a Zn2Cys6 binuc- 
lear cluster [25] - activates phol^ expression. Cskl - a 
CDK-activating kinase-activating kinase (CAKAK) [26] 
- represses phol^ expression in high-Pi conditions. Epista- 
sis analysis indicates that Pho7 acts downstream of Cskl. 

In this study, we explore how these factors affect tran- 
scriptional output by characterizing the PHO transcrip- 
tional response in S, pombe. We analyze this response as 
a function of phosphate, Pho7, and Cskl availability 
using DNA microarrays. We delineate a core PHO tran- 
scriptional response comprising the genes phol^, 
SPBC8E4.01C (an S. cerevisiae PH084 ortholog), and 
SPEC 127 1.09 (an S. cerevisiae GITl ortholog), whose in- 
duction in response to phosphate starvation is conserved 
between S. cerevisiae and S. pombe. Interestingly, while 
these three genes share a functionally analogous regula- 
tory pathway (i.e. activation through a transcription fac- 
tor that is normally repressed by a kinase) we find that 
the mechanism for regulation differs widely between 



species. Our analysis of the Pho7-regulated transcrip- 
tional output - coupled with a global profile of Pho7 
binding to promoters of stress responsive genes - leads 
us to the conclusion that, unlike Pho4, Pho7 plays a role 
in multiple stress response pathways. We conclude that 
while there is a core PHO transcriptional response 
shared between these two ascomycetes, the systems logic 
and specialization of PHO components varies widely. 

Results and discussion 

pho7^ and csk1^ regulate a core subset of the PHO 
response in S. Pombe 

The kinetics and maximal output of transcription vary 
widely between different environmental stress response 
pathways [27-29]. To outline the specific PHO response 
in S, pombe for subsequent analysis, and to avoid indir- 
ect activation of non-phosphate starvation regulated 
genes, we performed a single time-dependent, genome- 
wide expression analysis of wild-type S, pombe cells in 
medium lacking inorganic phosphate (no-Pi conditions, 
see Additional file 1). 

The starvation time-course revealed two distinct 
responses to phosphate starvation. The rapid response 
contained 63 genes that exhibited an increase in expres- 
sion at 120 minutes post-starvation (red lines in 
Additional file 1, genes listed in Additional file 2; see 
Methods for gene selection criteria). This class contains 
the secreted acid-phosphatase, phol^ (orthologous to 
ScPHOS), SPBC8E4.01C (orthologous to ScPH084), and 
SPBC1271.09 (orthologous to ScGITl), As phol^ and 
SPBC8E4.01C induction has been previously observed in 
response to Pi starvation [24,30], we believe that this set 
accurately reflects the genes that respond rapidly to 
changes in external Pi. In contrast, the slower response 
(86 additional genes, blue lines in Additional file 1, genes 
listed in Additional file 2; see Methods for gene selection 
criteria) was significantly enriched for genes previously 
implicated in a generalized stress response [27]. We fo- 
cused our attention on the fast responding genes to 
avoid indirect effects caused by persistent stress in cells. 

Previous work indicated that phoT^ and cskl^ are im- 
portant regulators of phol^ expression [24]; we expected 
that they would also play a significant role in regulating 
additional components of the PHO response. To test 
our hypothesis, we probed the transcriptional profiles of 
wild-type, pho7A, cskl A, and pho7AcsklA strains in 
high-Pi and no-Pi conditions at 120 minutes post- 
starvation using DNA hybridization microarrays. 

We identified 22 genes induced in response to Pi limita- 
tion (Figure lA, first column, applying a p-value cutoff of 
0.10; see Methods for gene selection criteria). If expression 
of these genes is dependent on the activity of phoT^, then 
their transcript abundance should decrease in a pho7A 
strain when compared to a wild-type background (both 
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Figure 1 Pho7 and Cski regulate a core set of genes induced 
in response to piiospiiate limitation. (A) Heat map showing the 
fold change (log2 scale) of Pi-starvation induced genes in a 
wild-type background (first column), a pho7A strain in high- and 
no-Pi conditions vs. a phoT^ strain (second and third), a csklA strain 
in high- and no-Pi conditions vs. a cskl^ strain (fourth and fifth), a 
pho7/^cskl /S. strain in high- and no-Pi conditions vs. a phoT^cskl^ 
strain (sixth and seventh), and a cs/c/A strain in high- and no-Pi 
conditions vs. a pho7AcsklA strain (eighth and ninth) as measured 
by microarray analysis. Thresholds used for gene selection are 
described in the materials and methods. For a full list of Pi, phoT^, 
and cskl^ dependent genes see Additional file 3. (B) Pie chart 
showing the percentage of genes induced in response to Pi 
limitation with Pi- and Pho4-dependent orthologs in 5. cerevisioe [7]. 



starved of Pi). Induction of 31.8% (7/22) of these genes 
during Pi starvation is phoT^ dependent (Figure lA, third 
column). If the repressor {cskl^) prevents this induction 
by phoT in high-Pi conditions, then transcript abundance 
in a csklls. strain compared to the cskV background in 
high- Pi conditions should increase. Three genes display 
this response (Figure lA, fourth column). A complete list- 
ing of all genes regulated by Pi, phoT^, and/or cskl^, along 
with their orthologs in S, cerevisiae can be found in Add- 
itional file 3. Finally, comparisons of the phoVb^, csklb^, and 
pho7^cskl^ double deletion strains confirm the previously 
described epistatic relationship between phoT and cskl^ (i. 
e. a loss of phoT in a cskl^ switches constitutive 



expression to uninducible expression) (Figure lA, last two 
columns). 

When phosphate is limiting, S, cerevisiae Pho4, along 
with Pho2, induces the transcription of genes required for 
phosphate acquisition [7]. The orthologs of Pho4, Pho2, 
Pho81, and PhoSO are not found in S, pombe or involved in 
the PHO response, raising the question: does a functionally 
analogous signaling pathway involving phoT and cskl^ 
regulate the PHO transcriptional response in S, pombel 
Comparing the Pho4, Pi-starvation induced genes in S, cer- 
evisiae [7] with the phoT , Pi-starvation induced genes in 
S. pombe reveals an overlap of only three orthologs 
(Figure IB). For phor, SPBC8E4.01c, and SPBC1271.09 a 
similar system of transcription factor activation, with re- 
pression by a kinase in high-Pi conditions, occurs. Unlike 
Pho85-Pho80 regulation of Pho4 in S. cerevisiae, most of 
the /^/zoT^-mediated response is independent of cskl^ regu- 
lation. Further, a large segment (82%) of the observed 
S. pombe PHO response is not conserved between S. pombe 
and S. cerevisiae. Therefore, the primary role for phoT and 
cskl^ in the PHO pathway is regulating the core set of 
phosphate harvesting and transport genes (phol^, 
SPBC8E4.01Q and SPBC1271.09). 

Pho7 is enriched at the PHO core promoters during Pi 
starvation 

Pho7 is classified as a putative transcription factor 
because it possesses a Zn2Cys6 binuclear cluster (ZC), a 
DNA binding domain for a number of transcription 
factors [31,32]. To test if Pho7 binds to j^/zoT^-regulated 
promoters, cells containing a functional, epitope 
tagged version of Pho7 (Pho7-TAP) were grown in 
high-Pi or no-Pi medium and purified DNA asso- 
ciated with Pho7 was processed via high-throughput 
sequencing (ChlP-Seq, see Methods). 

Surprisingly, there is widespread Pho7 binding even in 
high-Pi conditions (1676 peaks out of 4208 passed peak-ID 
thresholds; see Methods) (Additional file 4). During Pi star- 
vation 367 Pho7-bound sites exhibit an increase in Pho7 
enrichment. The highest levels of enrichment were 
observed in the promoters of Pho7-regulated genes identi- 
fied in the microarray analysis (for a complete list of the 
identified peaks, see Additional file 5). Further, there is a 
distinct overlap between genes whose expression levels are 
regulated by Pi and/or phoT and those whose promoters 
display enrichment in Pho7 binding. 13 of the 22 Pi 
dependent genes (p-value = 8.1e-8), 16 of the 29 phoT 
dependent genes (p-value = 9.5e-9), and 6 of the 7 genes 
regulated by both Pi and phoT (p-value = 1.3e-5) have 
promoters that are enriched for Pho7 binding in 
no-Pi medium (p-values were determined using a hyper- 
geometric test, see Methods for details). These results are 
very different from the global binding profile of Pho4 in S. 
cerevisiae. In that system, Pho4 is only recruited to the 
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promoters of PHO regulated genes during phosphate star- 
vation and, even then, only to relatively few locations (115) 
within the genome [7] . 

Pho7 binding is significantly enriched in the 
promoters of phol^, SPBC8E4.01c, SPBC1271.09, 
SPAC9E9.09C, SPBPB2B2.05, and SPCC330.06c in 
no-Pi conditions (Figure 2). In the case of phol^, 
SPBC8E4.01C, SPAC9E9.09C, and SPBPB2B2.05 there is 
enrichment of Pho7 even in high-Pi conditions (when 
compared to mock treated samples). This explains the pre- 
viously noted basal expression of phol^ in S, pombe [30,33] 
- it appears that Pho7 is bound and activating transcrip- 
tion even in the absence of stress (Figure lA, second 
column). Based on our microarray results, we note that 
increased Pho7 enrichment during Pi starvation is not due 
to an increase in phoT transcript abundance (0.05 log2 
fold-change in no-Pi vs. high-Pi conditions). Moreover, we 
observe no significant difference in Pho7-TAP protein 
levels during Pi starvation (Additional file 6). Together 
these observations suggest that the enrichment of Pho7 in 
the promoters of Pho7-regulated genes is attributable to 
increased affinity for the promoter and not due to an in- 
crease in Pho7 abundance. Further, the number of distinct 
Pho7 binding events varies between different promoters 
(Figure 2), and we do not observe a clear correlation 
between the number of binding events and transcriptional 
up-regulation. Our attempts to identify a DNA binding 



motif for Pho7 based solely on our ChlP-Seq data were 
unsuccessful (data not shown). 

Cski does not regulate Pho7 promoter occupancy 

Motivated by the following observations, we studied 
what effect loss of Cskl would have on Pho7 enrichment 
in high-Pi conditions: (1) Cskl represses the expression 
of the core PHO genes (phol^ , SPBC8E4.01c, and 
SPBC1271.09) in high-Pi medium (Figure IB, column 4); 
and (2) a decrease in Pi results in enrichment of Pho7 at 
the core PHO promoters (Figure 2). If binding of Pho7 
to the PHO promoters is necessary to drive increased 
transcriptional output, and Cskl represses transcription 
by preventing Pho7 binding, then a loss of Cskl in high- 
Pi conditions should result in an increase in Pho7 bind- 
ing (mimicking the no- Pi binding profile observed for 
Pho7 in wild- type cells). Pho7 binding at the phol^ pro- 
moter in a C5/dA background is mildly increased com- 
pared to a cskl^ background (Figure 3). This increase in 
binding is less than the observed increase in Pho7 bind- 
ing in cskl^ cells grown in no- vs. high-Pi conditions 
(1.2-fold versus 4-fold), suggesting that Cskl is not the 
major regulator of Pho7 binding at the phol^ promoter. 

We then examined the global effect of Cskl loss on 
the binding profile of Pho7 using ChlP-Seq with cskl A 
cells grown in high-Pi conditions. Unlike the enrichment 
during Pi starvation, deletion of Cskl does not result in 
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Figure 2 Pho7 is recruited to piiospiiate-regulated promoters during Pi starvation. Cells lacking the TAP epitope tag (Mock) or containing 
an epitope fused version of Pho7 (Pho7-TAP) were assayed for their enrichment of Pho7 at the phol^ promoter via ChlP-Seq. Cells were either 
incubated in high-Pi conditions (blue) or starved (red) for 2 hours before cross-linking. Reads were normalized to total counts for each 
chromosome. Magnified ChlP-Seq profiles for phol^, SPBC8E4.01c, SPBC1271.09, SPAC9E9.09c, SPBPB2B.05, and SPCC330.06c are shown. Mock 
samples lacking the TAP epitope (black) were treated identically to experimental samples. Coordinates are relative to the initial ATG, with gene 
products identified by black bars. 
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Figure 3 Cski does not regulate Pho7 Promoter enrichment. ChlP-Seq profiles for phol^, SPBC8E4.01c, SPBC1271.09, SPAC9E9.09c, 
SPBPB2B.05, and SPCC330.06c are shown. Cells from either a cskl^ (blue) or csklA (cyan) background were incubated in high-Pi conditions for 2 
hours. Mock samples lacking the TAP epitope (black) were treated identically to experimental samples. Coordinates are relative to the initial ATG, 
with gene products identified by black bars. Reads were normalized as previously described. 



a global increase in Pho7 binding in high-Pi conditions 
(Additional file 7). At the core PHO responsive genes we 
observe either no change (SPBC8E4.01c) or a slight in- 
crease in Pho7 binding in the csklA strain (phol^ and 
SPBC1271.09), which is still well below the enrichment 
seen during Pi starvation (Figure 3). As we observed 
during Pi starvation, the loss of Cskl does not influence 
either phoT transcript abundance (-0.12 log2 fold- 
change, [c5/dA vs. cskl^\ -nPi) or Pho7-TAP protein 
levels (Additional file 6). We draw two conclusions from 
these data: (1) the level of Pho7 bound in high-Pi condi- 
tions would be sufficient to induce high levels of tran- 
scription if not for the repressive action of Cskl; and (2) 
Cskl does not repress Pho7 activity by preventing Pho7 
from binding to the promoters of responsive genes. 

A Pho7 upstream activating sequence (UAS) and an 
independent Pi sensing module control phol^ expression 

Based on our ChlP-Seq results, we know Pho7 binds 
between nucleotides -280 and -180 in the phol^ promoter. 
To determine whether the sequences in this region are ne- 
cessary and/or sufficient for Pho7-dependent, Pi-starvation 
induced expression, we utilized an in vivo strategy for con- 
firming Pho7-promoter interactions using exogenous ex- 
pression plasmids. Briefly, differing lengths of the phol^ 
promoter driving the expression of yellow fluorescent pro- 
tein (yfp^) were constructed in a vector and transformed 
into phoT, pho7Ay and cskl A backgrounds (Figure 4A). 
yjp"" expression was measured using FACS with mean YFP 



intensity serving as a proxy for promoter activity (see 
Methods). The 2 kb segment of the phol^ promoter acti- 
vates yfp"" expression during Pi-starvation - it exhibits a 
~ 7-fold increase in YFP intensity upon starvation 
(Figure 4B). This induction is dependent on Pho7, as it is 
abolished in a pho7A background (Figure 4B). Surpris- 
ingly, we find that loss of Cskl does not result in high- 
levels of yfp"" expression with the 2 kb reporter in cskl A 
cells grown in either high-Pi or no-Pi growth media. 
While the levels of expression are above those seen in 
phoT cells grown in high-Pi conditions and pho7A cells 
grown in any condition, they are significantly lower than 
the levels observed in a phoT background in no- Pi condi- 
tions (Figure 4B). Based on our evidence that trimming 
the phol^ promoter beyond 1 kb results in high level, 
constitutive expression of yfp^ in cskl A ceUs (Figure 4D, 
discussed below), our results with the 2 kb fragment in 
the cskl A background might be complicated by: (1) an 
additional repressor element located between 2 kb and 
280 bp in the phol^ promoter; and/or (2) differences be- 
tween the chromatin structure/promoter architecture of 
the endogenous phol^ locus and the exogenous 2 kb phol 
^P^'yJp^ vector that are influencing transcription in the 
cskl A cells. Therefore, we focused on the behavior of the 
shorter construct (280 bp phol^pr-yjjp^; Figure 4D) as a 
proxy for the interaction of Cskl at the region bound by 
Pho7-TAP in the phol^ promoter. 

If the phol^ promoter sequence bound by Pho7-TAP 
in the ChlP-Seq experiment is necessary for activation of 



Carter-O'Connell et al. BMC Genomics 2012, 13:697 
http://www.bionnedcentral.conn/1 471 -21 64/1 3/697 



Page 6 of 1 5 



I yfp+ 

2 kb pho1+pr-yfp+ 



UASA phoi yr-yfp+ 

\ yfp+ 



280 bp pho1 +pr-yfp+ 

\ yfp+ I 

180 bp plio1+pr-yfp+ 



6000 




]+Pi[ 



YFP Logio(A.U.) 



pho7+ pho7A csk1A 
2 kb pho1+pr-yfp^ 



8000 



^6000 
< 

S: 4000 



i 2000 



8000 



^6000 
< 

S: 4000 
> 



i 2000 



pho7+ pho7A 



E \JASA pho1+pr-yfp+ 



pho7^ pho7A csk1A pho7^ pho7A csk1A 
280 bp pho 7+pr-y^p+ 1 80 bp phol'^^r-yfp'^ 







9? / ^ 

' '/ ??^^/Pho7) 


/ ??,,/Plio7J 


■g Pi Sensing UAS ■■^^B 


1 pho1+ 1 


-2kb -280 bp -180 bp 


-2kb -280 bp -180 bp 



Figure 4 A Pho7 UAS in thie phol^ promoter is necessary expression. (A) Schematic of pho 1^ promoter variants used for expression analysis. 
2 l<b phol^pf-yfp'^: full promoter, 280 bp p/io/^pr-yfp^: promoter trimmed to only the Pho7-UAS, 180 bp pho/^pr-yfp^: promoter trimmed beyond 
the Pho7-UAS, UASA phol^pf-yfp^: full promoter lacking only the 20 bp under the Pho7 ChlP-Seq peak. See text for vector construction. (B) FACS 
histogram profile for the 2 kb phol^pr-yfp'^ construct transformed into a phoT^ strain. Cells were grown in either high-Pi (blue) or no-Pi (pink) 
media for 4 hours prior to fixation and counting. Bar graph depicting mean YFP intensity for the 2 kb phol^pr-yfp'^ construct transformed into a 
phoJ^, pho7A, and csklA strain and grown in either high-Pi (blue) or no-Pi (pink) media are given as the average of three biological replicates ± 
SE. (C) The UASA phol^p(-yfp^ construct was transformed into phoZ^ or pho7A strains and treated as in B. (D) The noted constructs were 
transformed into phoZ^, pho7A, or csklA backgrounds and treated as in B. (E) Model detailing the role of Cskl in repressing Pho7 function at the 
phol^ promoter. In high-Pi conditions Pho7 is not recruited to the UAS and Cskl interacts near the UAS to prevent full Pho7 activity. When the 
external concentration of Pi drops, additional Pho7 is recruited to the promoter and Cskl repression is relieved resulting in maximal induction. 
Proper regulation is controlled within a Pi sensing module located between bases -2000 and -280 in an unknown manner. 



phol^ transcription during Pi-starvation, then deletion 
of this region should result in a loss of y/p"" expression 
during Pi-starvation. To test this hypothesis, we gener- 
ated a construct in which the 20 bp centered under the 
Pho7-TAP ChlP-Seq signal were deleted (Figure 4A, 
UASA phol^pr-yfp^). In high-Pi growth conditions, the 
loss of the Pho7 bound region results in a slight increase 
in yfp^ expression (Figure 4C, compare blue columns to 
the blue columns in Figure 4B), and in Pi-starvation this 
UASA phol^pr-yfp^ construct no longer fully activates 
yfp^ expression (Figure 4C, compare column 2 to column 
2 in Figure 4B). The loss of phoT^ results in a further de- 
crease in expression from this UASA phol^pr-y^^ 
construct (Figure 4C, column 4). It is possible that Pho7 
recognizes additional segments of the promoter, though 



such contributions to activation in Pi-limiting conditions 
are modest. Together these results demonstrate that the 
Pho7-TAP bound promoter element is necessary for 
Pho7-dependent transcriptional activation during Pi-limi- 
tation. We have termed this region the Pho7 upstream ac- 
tivating sequence (UAS). 

To test whether the Pho7 UAS is sufficient for Pho7- 
dependent. Pi-limitation induced transcriptional activation 
we deleted all but the first 280 bp of the phol^ promoter 
and assayed in vivo yjjp^ expression (Figure 4A, 280 bp 
phol^pv-yfp^), yjp^ expression from the 280 bp phol^ 
pv-yfp'' construct is elevated in high- Pi conditions, and is 
only marginally activated during Pi-starvation (Figure 4D, 
phoT, column 1). Expression from this construct is 
reduced in a pho7l^ background and is unaffected by Pi 
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limitation - thus, expression in high-Pi conditions and the 
modest expression increase in Pi-Umiting conditions are 
dependent on Pho7 (Figure 4D, pho7A, column 2). In each 
background tested, the mean YFP intensities from the 280 
bp phol^pr-yjjp^ construct vary by less than 1.5 fold be- 
tween high-Pi and no- Pi conditions (Figure 4D, left panel). 
This is in contrast to both yfp^ expression from the 2 kb 
phol^pr-yjjp^ construct and endogenous expression of 
phol^ during Pi-starvation which exhibit > 10-fold induc- 
tion in Pi limitation. Thus, the Pho7 UAS is necessary but 
not sufficient for Pho7-dependent transcriptional activa- 
tion during Pi-starvation. 

Interestingly, the 280 bp phol^pv-yfp^ construct is cap- 
able of inducing full yfp"" expression in no-Pi conditions in 
a csklA background (Figure 4D, csklA, column 3). How- 
ever, the 280 bp phol^pv-yfp^ construct is not capable of 
relieving Cskl repression in Pi-starvation conditions - ex- 
pression is not induced in response to Pi limitation. Trim- 
ming beyond the Pho7 UAS results in transcriptionally 
inactive promoters in all backgrounds tested (Figure 4D, 
columns 4-6, 180 bp phol^pr-yjjp^). We conclude that 
there must be promoter elements present in the region 
between -2 kb and -280 bp in the phol^ promoter that act 
as a Pi-sensor: (1) preventing partial Pho7-dependent acti- 
vation in high- Pi conditions; and (2) are important for 
Cskl de-repression during Pi-starvation. 

Our FACS and ChlP-Seq results lead us to the follow- 
ing model for Pho7 and Cskl regulation at the phol^ 
promoter. In high-Pi conditions some Pho7 is bound to 
the UAS in the phol^ promoter. Pho7 in this state drives 
basal expression of phol^, Cskl, through an interaction 
with either Pho7 or elements near the UAS (directly or 
indirectly), prevents the full activation of phol^ expres- 
sion. The upstream Pi-sensor in the promoter ensures 
that Cskl remains repressive in these conditions through 
an as yet unspecified mechanism. During Pi starvation, 
Cskl repression is relieved and additional Pho7 is 
recruited to the phol^ promoter, driving maximal 
expression (Figure 4E). Investigating the promoter ele- 
ments and transcription factors that comprise the 
Pi-sensor - as well as the use of this promoter structure 
at additional Pho7-dependent and Pi-starvation indu- 
cible promoters - is an exciting area for future research. 

Pho7 regulates gene expression in response to multiple 
stress conditions 

During our expression analysis we noticed a set of genes 
with decreased expression in high phosphate conditions 
in a pho7A background that are not induced during phos- 
phate starvation (Additional file 3). Additionally, Pho7 is 
bound to the promoters of a number of these genes in the 
ChlP-Seq analysis (Additional file 8). These observations 
raise the following question: is Pho7 dedicated solely to 



the phosphate starvation pathway, like Pho4, or does it 
play a broader role in the stress response? 

To answer this question, we identified the genes from 
our microarray analysis that display an increase in expres- 
sion between the wild-type and pho7A backgrounds (in 
either +Pi or -Pi conditions; thresholds described in mate- 
rials and methods) and asked whether their promoters 
contained a significant peak of bound Pho7 within 800 bp 
of the start codon - there are 63 genes that meet these 
requirements. This gene set was then processed through 
the Gene Ontology Tools: Term Enrichment algorithm 
[34], allowing us to look for biological functions over- 
represented in our gene set. For a GO term to be consid- 
ered enriched, we required at least three unique gene pro- 
ducts be included in the term family and the enrichment 
must have a p-value < 0.01. 32 of the original 63 genes 
meet these requirements and are classified by the highest 
parent term available (Figure 5A). 15 of these 32 genes 
are involved in transmembrane transport, including 3 
involved in phosphate ion transport. Characterizing the 
transmembrane category in further depth we find genes 
involved in iron, copper, and zinc transport (Figure 5B). A 
complete listing of GO Terms and the enriched gene set 
can be found in Additional file 9. 

From the 32 identified genes, we utilized gpdl^, hxk2^y 
fiol^y and ctr4^ expression as proxies for Pho7-mediated 
transcriptional induction in various stress conditions. 
Gpdl is a glycerol-3-phosphate dehydrogenase that 
synthesizes glycerol and is essential for survival during 
osmotic stress [35]. During osmotic stress, glycerol pools 
increase, protecting the cell. Hxk2 is a hexokinase that 
plays a role in regulating alternative carbon utilization 
when glucose sources are limited [36]. It is maximally 
induced in response to a switch from glucose to glycerol 
as a carbon source. Fiol, in conjunction with Fipl, com- 
prises the oxidase-permease iron transport system 
responsible for harvesting iron in depleted conditions 
[37]. During iron repletion fiol^ is repressed by the 
activity of Fepl, an iron sensing transcription factor [38]. 
During iron starvation fiol^ is de-repressed and induced 
~70-fold. Finally, Ctr4 is a high-affinity copper trans- 
porter that is induced in copper depletion conditions by 
Cufl, a copper sensing transcription factor [39]. We 
designed RT-qPCR primer sets for each of these genes 
and measured their expression as a function of osmotic, 
iron, copper, and carbon utilization stress in both phoT" 
and pho7A backgrounds (see Methods). 

As previously demonstrated [24], loss of Pho7 com- 
pletely abrogates induction of phol^ in no-Pi medium 
(Figure 5C, first panel). For each of the additional stresses 
tested, the loss of Pho7 causes a significant decrease 
(p-value < 0.05) in the maximal induction of the target 
gene (Figure 5C). The Pho7 dependence of these genes 
varies, with some ifiol^y ctr4^) showing a relatively minor 
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Figure 5 pho7^ functions in additional stress responses. (A) Genes with a Pi- and p/ioAcomponent based on microarray analysis were cross- 
referenced with a list of genes with phoT^ enrichment within 800 bp of the initial ATG (ChlP-Seq). The pie chart shows the GO Term enrichment 
for the identified genes as compared to a background model of all 5. pombe transcripts. The boxed inset highlights a number of genes that were 
found in the small molecule catabolic process class, with those in bold serving as the initial targets for expression analysis. P-values for the 
enrichment against the background model are given. (B) The transmembrane transport category from A was expanded to identify specific ions 
whose transport is regulated by phoT^. Bold candidate genes served as initial target for expression analysis. (C) phoT^ cells were incubated in 
either normal (white) or stress (light gray) conditions, and the expression of the appropriate gene was measured by RT-qPCR. pho7A cells were 
also stressed (dark grey) to identify the level of p/io^-dependent regulation. Expression was normalized to octl^. Shown is the average of three 
biological replicates ± SE. (D) Heat map displaying the fold-change (log2 scale) for genes induced during stress in a p/io/^-dependent manner. 
The top track in each sub-column compares expression in replete versus stress conditions for wild-type cells; the bottom track compares phoZ^ 
to pho7A cells in the displayed stress (Pi: phosphate starvation, Fe: iron starvation, Cu: copper starvation, NaCI: osmotic shift, GE: carbon switch). 
Boxed regions identify genes that are induced in response to stress and require phoT^ for their stress response. Pie charts depict the percentage 
of genes responding to stress in a phoT^ independent manner (light grey), repressed by the loss of phoT^ but unaffected by stress (dark grey), 
and dependent on both (red). 



Pho7 component, while others (gpdl^y hxk2^) appear fully 
dependent on Pho7 to reach a peak level of induction in 
stress. None of the Pho7-regulated responses were as 
dramatic as that observed for Pi starvation, which may in- 
dicate that Pho7 plays a more subtle role in coordinating 
expression at Pi-independent loci. Given this subtle re- 
sponse, the only recent availability of a deletion collection 
for S, pombe [40], and the fact that Pho7 function was 



unverified until recently, it is not surprising that the more 
general role for Pho7 has not been previously observed. 

It remains possible that the stress response effects we 
observe are artifacts limited to only the few genes we 
studied using RT-qPCR. Using microarray analysis with 
RNA collected from pho7^ or pho7A cells grown in non- 
stress and stress conditions, we examined the stress 
response mediated by phoT in the above conditions. 
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We find that of 274 genes induced in the stress conditions 
studied, 44 genes are phoT dependent (7 genes are 
induced in multiple stress conditions) (Figure 5D and 
Additional file 10). phoT is responsible for coordinating 
between 11.7-23.5% of the total stress response in each 
condition (for comparison, phoT coordinates 21.7% of 
the Pi starvation response using these thresholds). In each 
stress we find enrichment of distinct GO terms. The set of 
iron-responsive, /^/zoT^-regulated genes is significantly 
enriched for the biological process of iron assimilation 
(Jiol", fipr, str3^, sib2^, p-value = 5.2e-09), the set of 
copper- responsive, jE^/zoT^-regulated genes is enriched for 
copper ion transport {ctr4^ and ctrS^, p-value = 5.6e-05), 
and the osmotic shock-responsive, /^/zoT^-regulated gene 
set is enriched for metal ion transport {zrtl^ and ctr4^, 
p-value = 7.8e-04). With the carbon switch, /^-regu- 
lated response we see an enrichment of genes responsible 
for small molecule metabolism as well as conjugation. The 
reasoning for the conjugation process enrichment or why 
phoT would be involved is unclear. Overall, the biological 
processes are generally linked with transmembrane trans- 
port, suggesting that phoT is responsible for coordinating 
the correct transport of nutrients required for each stress 
response. Unlike the system in S. cerevisiae, where the cen- 
tral Pi-starvation regulator is tightly linked with the PHO 
response, the phoT based system in S. pombe functions 
differentially in a number of stress response networks. 

Conclusions 

In this study we have defined and characterized the gene 
regulatory network in S, pombe responsible for coordin- 
ating the response to inorganic phosphate starvation. 
There are two distinct temporal responses in the PHO 
pathway in S, pombe: a fast response concerned with im- 
mediately harvesting inorganic phosphate from the en- 
vironment and transporting it into the cell, and a slower 
one associated with a general stress response. Within 
the fast response we define a core PHO regulon com- 
prised of the phol^, SPBC8E4.01C, and SPBC1271.09 
genes whose induction in response to phosphate starva- 
tion, and regulatory behavior, has been conserved be- 
tween S, pombe and S, cerevisiae. 

Unlike the PHO response in S. cerevisiae, however, the 
positive regulator in S. pombe is bound to the phol^ and 
SPBC8E4.01C promoters irrespective of external phos- 
phate availability. Using our Pho7-TAP ChlP-Seq data- 
set, we identified a single Pho7 binding region in the 
phol^ promoter. This region - located between nucleo- 
tides -265 and -245 - serves as a Pho7-dependent UAS, 
which is necessary for transcriptional activation of phol^. 
The interaction of Pho7 with the UAS leads to a basal 
expression of the secreted acid phosphatase in high-Pi 
conditions. The 1.8 kb region preceding the UAS in the 
phol^ promoter is required for Pi sensing, coordinating 



the activation by Pho7 and repression by Cskl based on 
Pi availability. Cskl prevents full activation of Pho7 during 
phosphate replete conditions and repression is maintained 
even in the minimal UAS construct (280 bp phol^pr-yjp^). 
During phosphate starvation this inhibition is relieved 
(through an unknown mechanism) and additional Pho7 is 
recruited to a number of sites throughout the genome, 
causing further induction of /^/zoT^-dependent genes (see 
Figure 6 for a comparison of the S. cerevisiae and S, pombe 
PHO pathways). In previous work, Cskl was shown to 
regulate transcription by activating the positive transcrip- 
tion elongation factor b (P-TEFb) ortholog, Cdl<9 [41]. 
Cdl<9 coordinates transcript elongation and processing, 
and its full activation by Cskl leads to an increase in CTD 
kinase activity [42]. How this generally positive regulatory 
network is switched to an inhibitory role in the PHO sys- 
tem remains an open question. 

We were also surprised to find that Pho7 was bound 
throughout the genome in both high-Pi and no-Pi condi- 
tions. We had thought based on previous evidence that 
Pho7, like Pho4, would be specific to the PHO response. 
Instead we demonstrate that Pho7 binds within the pro- 
moters of additional stress responsive genes and plays a 
role in iron, copper, osmotic, and alternative carbon 
utilization stress. Each stressor elicits a different pho7 
^-dependent transcriptional response, though it appears 
that the main regulatory role of phoT is coordinating 
stress-specific transmembrane transport. There must 
exist some mechanism to either direct Pho7 to the 
proper location for inducing the correct genes or acti- 
vate Pho7 at only the appropriate locations (or some 
mixture of both). In S. cerevisiae, the osmotic, oxidative, 
and glucose limitation stress responses are mediated by 
the transcription factor Msn2 [43]. In normal conditions, 
Msn2 is phosphorylated and its entry into the nucleus is 
limited [44]. Different stresses elicit distinct dynamics of 
nuclear transport, leading to different transcriptional 
outputs [29]. Given that Pho7 is bound to the genome 
constitutively, we do not expect that nuclear exclusion 
will play as large a role as it does with Msn2 regulation, 
but it remains possible that differential post-translational 
modifications are responsible for this combinatorial acti- 
vation by Pho7. Pho7 may be playing a more passive role 
in regulation, with additional factors determining Pho7 
genomic localization. 

Nonetheless, we have demonstrated that within the 
evolutionary parallel signal transduction networks that 
comprise the PHO pathway there exists a core PHO 
transcriptional regulon. The specific mechanisms 
involved in regulating the PHO response in S. cerevisiae 
and S. pombe show remarkable flexibility. An interesting 
area for future research centers on the environmental 
factors that contributed to the development of these two 
parallel networks. Why is the PHO response in S. pombe 
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Figure 6 A comparison of the S. cerevisiae and S.pombe PHO signaling pathways. (A) A schematic depicting tine PHO signaling patliway in 
5. cerevisiae. In conditions where Pi is plentiful (high Pi, left) Pho4 is multiply phosphorylated by the Pho85-Pho80 complex and is excluded from 
the nucleus. Depletion of inorganic phosphate (no Pi, right) leads to the IPy-PhoSl -dependent inhibition of Pho85-Pho80. Relief of Pho85-Pho80 
mediated repression allows unphosphorylated Pho4 to enter the nucleus, interact with Pho2 at PHO-specific genes, and induce transcription. 
(B) A schematic depicting the PHO signaling pathway in 5. pombe. In high Pi conditions (left) some Pho7 is bound to the promoters of Pi- and 
pho^-dependent genes. Bound Pho7 is prevented from reaching maximal activation by the repressive effects of Cskl. Loss of external phosphate 
(right) results in the recruitment of additional Pho7 to the promoter and the relief of Cskl -mediated repression. The Pho7-dependent genes 
identified in this study that are regulated in this manner are listed (5. cerevisiae orthologs in parenthesis). 



under the control of a general stress transcription factor, 
Pho7, while S. cerevisiae has developed the phosphate 
starvation specific pathway for Pho4 activation? What 
are the environmental pressures that favor a "leaky" 
response in S, pombe and a tightly controlled one in S, 
cerevisiae! Broadly speaking, our study provides a frame- 
work for determining the fundamental requirements for 
regulating phosphate homeostasis in Ascomycota and 
the specific points in the signal transduction pathway 
that can be altered as conditions merit. 

Methods 

Growth conditions and strains 

S, pombe cells were maintained in previously described 
YES or EMM media [45]. The yeast strains used were: 
DPI (972 /z), DP18 {ura4-D18 ade^ leu" h^), DP81 
(pho7AI<ANMX6 972 /z ), DP94 (Pho7-TAP/04A/MX6 



972 /z), DP106 {csklANATMX6 972 /z), DP109 
{csklM<ANMX6 ura4-D18 ade" leu" DPI 11 

(pho7M<ANMX6 ura4-D18 ade" leW h^), DPI 13 
(pho7AKANMX6 csklANATMX6 972 /z ), DP114 (Pho7- 
TA]?KANMX6 csklANATMX6 ura4' 972 /z ), and DP115 
(Pho7-TAP/047VMX6 csklANATMX6 ura4- 972 h^). The 
functionality of the Pho7-TAP allele was confirmed by 
both liquid phosphatase assay and RT-qPCR analysis 
and it behaves as pho7^. To tag Pho7 and delete cskl^ 
we utilized a PGR fragment containing the marker 
of interest flanked by homologous regions for the spe- 
cific gene target. Cells were transformed with lithium 
acetate and polyethylene glycol 8000 [45]. Primers used 
for deletion or tagging are found in Additional file 11. 
To provide consistency with previously published re- 
sults for inorganic phosphate starvation, all starvation 
experiments were conducted with cells incubated in a 
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90%SD-10%EMM media, which has been previously 
described [24]. 

Microarray analysis and data processing 

Strains were grown in 90%SD-10%EMM medium con- 
taining 10 mM KH2PO4 (high-Pi) at 30°C until they 
reached early-log phase (OD6oo=0-15-0.25). Cells were 
collected via filtration, washed twice, transferred to fresh 
media lacking Pi (no-Pi), and grown at 30°C for up to 4 
hours. Immediately prior to starvation (t=0), 20 mL of 
cells were added to 30 mL of methanol kept at -65°C to 
prevent further transcription or RNA degradation. At 
30, 60, 120, and 240 minutes post-starvation this process 
was repeated. Cells were left in methanol for 10 minutes, 
pelleted, washed quickly in autoclaved water, re- 
suspended in 750 uL of RNAlater (Ambion), and snap- 
frozen in liquid nitrogen. RNA was extracted using the 
RNeasy Mini kit (Qiagen). cDNA was generated in a re- 
verse transcriptase reaction using 10 [ig total RNA with 
a 1:1 mixture of oligo-dT and random hexamer primers 
(Operon) and a 2:3 ratio of amino-allyl-dUTP:dTTP 
(Sigma). Superscript II RT (Invitrogen) was added and 
the reaction mixture was incubated at 42°C for 2.5 
hours. cDNA was purified using a PCR purification kit 
(Qiagen) after completing hydrolysis of remaining RNA. 
An equal amount of cDNA from each time point was 
pooled to provide the reference sample. Purified cDNA 
samples were labeled using N-hydroxyl succimamide 
esters of either Cy3 or Cy5 dyes (GE Biosciences). 300 
ng of the Cy3 (each individual time point) and 300 ng of 
the Cy5 (pooled reference) labeled sample was competi- 
tively hybridized to custom Agilent 8x1 5K S. pombe two- 
color expression microarrays (GEO Platform:GPL15827) in 
2xGEx Hybridization Buffer (Hi-RPM) (Agilent) for 17 
hours at 60°C. Microarrays were washed and immediately 
scanned using an Axon 4000B scanner [46]. The mean in- 
tensity of each spot in the Cy3/Cy5 channels was extracted 
using the GenePix 5.1 software, followed by lowess and 
quantile normalization performed with the MATLAB 
bioinformatics toolbox. Expression ratios for each time 
point, X, were normalized to t=0 (Log2[Cy3t=x/Cy5pooi] - 
Log2[Cy3t=o/Cy5pooi]) and thresholds for induced genes 
were set at > 2a + median log2 fold change for each time 
point (1.00 log2 fold change at 120 minutes, 1.24 log2 fold 
change at 240 minutes). Genes above threshold at both 120 
minutes and 240 minutes post-starvation were classified as 
the rapid response. Genes above only the 240-minute 
threshold were classified as the slow response. The starva- 
tion time course was not repeated. Results for all of the 
microarray experiments conducted in this study are avai- 
lable through NCBI-GEO [GEO:GSE39478]. 

To determine the extent of phoT and cskl^ regulation 
within the PHO response we grew the relevant strains 
(DPI, DP81, DP106, DPI 13) as described above, with 



the exception that cells were split into either high- Pi or 
no-Pi media and grown for 2 hours prior to RNA collec- 
tion. Two independent biological replicates were per- 
formed for each of the conditions tested except for the 
phoTcsklMphoVt^csklb^ comparison in no-Pi media. 
For each of these arrays the two probes used to detect 
each ORE were averaged and treated as single data 
points with p-values determined using a students t-test 
with a one-tailed distribution against the null hypothesis 
in the MATLAB software. Thresholds were set at > 1.8 
log2 fold-change to facilitate comparison with the previ- 
ously characterized S. cerevisiae data set [7]. Genes pas- 
sing the induction threshold also had to pass a p-value 
threshold of < 0.10. Clustering analysis was completed 
using k-means clustering in the Cluster 3.0 program [47] 
after empirically determining the optimal number of 
clusters using the MATLAB bioinformatics toolbox. 

Chromatin immunoprecipitation of Pho7-TAP with 
high-throughput sequencing (ChlP-Seq) 

ChlP-Seq was performed on the DPI, DP94, and DPI 15 
strains as previously described [7,48-50]. Cells were 
grown to early log-phase (ODeoo^O-lS) in high-Pi media 
at 30°C and split into either 200 mL of high-Pi or no-Pi 
media and grown for 2 hours. Formaldehyde (Sigma) 
was added to a final concentration of 1% (v/v) to cross- 
link chromatin, and the reaction was allowed to proceed 
for 15 minutes. Glycine (Sigma) was then added to a 
final concentration of 125 mM and incubated for 5 min- 
utes to quench cross-linking. Cells were lysed by bead 
beating (6x2 min on, 2 min off) and chromatin was 
sheared to 300-600 bp fragments using a Misonix 
Sonicator 3000. Immunoprecipitation was performed 
with 100 uL of Protein G Dynabeads (Invitrogen) 
coupled to 4 uL of anti-Protein A antibody (Sigma). Pro- 
tein concentrations were measured using a Bradford 
Assay (Bio-Rad). Following the generation of ChIP lysate 
three aliquots of 650 (ig soluble protein were subject to 
immunoprecipitation and pooled just prior to elution 
from the beads. Samples were processed following the 
Illumina HTS guidelines with libraries of 200-300 bp 
selected via 2% agarose DNA gels. Libraries were ampli- 
fied by PCR and purity was determined using an Agilent 
High-Sensitivity DNA kit on an Agilent Bioanalyzer. 
Libraries were sequenced on an Illumina HighSeq 2000 
and 50bp reads were aligned to the S, pombe 972 h' ge- 
nome using ELAND. We obtained between 16 million 
and 54 million reads on average from our samples. 
Uniquely aligned reads were extended 80 bp from the 
read start site to cover the average length of insert as 
determined by the Agilent Bioanalyzer. Results for all of 
the ChlP-Seq experiments conducted in this study are 
available through NCBI-GEO [GEO:GSE39498]. 
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To determine which of our enriched regions were ac- 
tually attributable to a Pho7-TAP binding event we used 
a modified method from [48]. For each condition ana- 
lyzed we set a lower threshold for peak discovery equal 
to the genomic average of reads per base. We set the 
upper threshold equal to the highest observed read 
count within the given sample. Using 380 equal incre- 
ments between these two thresholds we defined peaks 
that were larger than 100 nucleotides and separated by 
at least 20 nucleotides. Peaks were compiled at the high- 
est threshold at which they met those standards and full 
peaks were required to be at least 150 nucleotides dis- 
tant from the nearest neighbor. Statistical analysis com- 
paring sample enrichment to mock enrichment was 
performed in MATLAB using previously described 
methods [48]. Peaks used for subsequent analysis had a 
> 2-fold enrichment over the genome average and a p- 
value < 0.005 compared to the mock sample. 

To determine the likelihood that the genes determined 
from microarray analysis to be regulated by Pi and/or 
pho7 would also have promoters that contain at least 
one Pho7 binding site we utilized a hypergeometric test. 
In this case the p-value is given by: 



where N is the total number of genes probed by the 
microarray (5046), M is the number of peaks within 
-800 and 0 bp of any start codon in no-Pi conditions 
(570), n is the full set of regulated genes (22 for Pi, 29 
for phoT^y and 7 for both), and k is the set of regulated 
genes with at least one Pho7 peak in the promoter (13 
for Pi, 16 for phoT , and 6 for both). 

Western blotting with Pho7-TAP 

DPI, DP94, or DPI 14 Cells were grown to log-phase in 
high- Pi media at 30° C. Following collection of a high- Pi 
sample, cells were washed three times in no- Pi media, 
transferred to no-Pi media, and grown for either 60 or 120 
minutes. Cells were lysed by beadbeating in urea lysis buf- 
fer (20 mM Tris-HCl, pH 8.0, 50 mM Na2HP04, 8M urea, 
and 1 mM PMSF). Total protein was quantified on a 
Nanodrop 2000 (Thermo Scientific) using a BSA standard 
(Bio-Rad). Equal amounts of total protein for each sample 
(30 ug) were subjected to separation by SDS-PAGE and 
transferred to nitrocellulose. Immunoblotting was per- 
formed with rabbit IgG (1:1000, Jackson ImmunoRe- 
search) followed by incubation with goat anti-rabbit-HRP 
(1:5000, Thermo Scientific). Blots were developed using 
the SuperSignal West Femto Chemiluminescent Substrate 



(Thermo Scientific) and analyzed on an Alpha Innotech 
Gel Imagining System. 

pho^^ promoter deletion analysis 

Segments of the phol^ promoter were amplified using 
PCR and cloned into a y/p"" plasmid using homologous 
recombination [51] containing the selectable ura^ mar- 
ker, creating a phol^pr-yfp^ fusion. Plasmids were trans- 
formed into DP 18, DP 109, or DPI 11 backgrounds using 
lithium acetate and polyethylene glycol 8000. Cells were 
selected based on their ability to grow in EMM-ura 
media. Cells containing the various plasmids were grown 
to early log-phase (ODeoo^O-lS) in high-Pi media (lack- 
ing uracil) at 30°C, collected, washed twice in sterile 
water, and split into either high-Pi or no-Pi media (both 
lacking uracil). Cells were grown for 4 hours at 30°C and 
100 uL of 10% buffered formalin (Sigma) was added to 
900 uL culture. Fixation proceeded for 5 minutes at 
room temperature prior to washing: once with O.IM 
potassium phosphate buffer (pH 8.5) and once with 
1.2M sorbitol in O.IM potassium phosphate buffer (pH 
8.5). Cells were resuspended in 1.2M sorbitol, O.IM 
potassium phosphate buffer (pH 8.5) and incubated 
overnight at 4°C. 

FACS counting with each sample was performed using 
a LSR II Analyzer (BD Biosciences). DP18, DP109, and 
DPI 11 cells lacking the yfp^ expression system were 
used to normalize forward-scatter, side-scatter, and auto- 
fluorescence for each experiment. 50,000 cells were 
counted for each experimental condition tested; cells 
with forward- and side-scatter values between 50,000- 
150,000 and YFP expression > mean autofluorescence 
were subject to further analysis. Three biological repli- 
cates were performed and the average YFP intensity for 
the replicates is reported ± SE. 

pho7^ regulation in additional stress response pathways 

For each individual stress response, initial cultures of 
DPI and DP81 were grown in 90%SD-10%EMM media 
containing 10 mM KH2PO4 to early-log phase at 30°C. 
Cells were collected, washed twice with autoclaved 
water, and split into the following conditions (all modifi- 
cations of the high-Pi media): +Pi - 10 mM KH2PO4, 
-Pi - 0 mM KH2PO4, +Fe - 100 uM Fe(III)Cl3, -Fe - 
250 uM 2-2'-dipyridine (DIP) (Sigma), +Cu - 100 uM 
Cu(II)S04, -Cu - 100 uM bathocuproine disulphonate 
(BCS) (Sigma), Osmotic Shift - 1.2M NaCl instead of 
O.IM NaCl, and Carbon Switch - 2% glycerol/1% etha- 
nol (GE) instead of 2% glucose (G). Cells were grown for 
2 hours and harvested as described above. Recovered 
RNA was converted into cDNA using the iScript cDNA 
synthesis kit (Bio-Rad) and subjected to RT-qPCR. Amp- 
lification of the gpdl^, fiol^, ctr4^, hxk2^, and phol^ 
transcripts were measured for the three independent 
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replicates and transcript abundance was normalized to 
actl^. Shown is the average ± SE. Primers used in the 
RT-qPCR analysis can be found in Additional file 11. 

Extracted RNA was also subjected to microarray ana- 
lysis as detailed above. Expression from phoT^ cells in 
replete conditions was compared to that in stress condi- 
tions for each individual stressor to determine the base 
set of genes that respond in each given stress. The 
dependence upon phoT was determined by comparing 
the levels of induction in phoT^ cells in stress to induc- 
tion in pho7A cells in stress. Based on previous reports 
[27] the osmotic shift conditions were assayed 20 min- 
utes post-shift to provide a more accurate measure of 
genes directly induced by osmotic pressure. Extraction 
of Cy3-Cy5 fluorescence intensity was performed using 
the GenePix 5.1 software and normalization was com- 
pleted using the MATLAB bioinformatics toolbox. At 
least two independent biological replicates were per- 
formed for each of the conditions tested. For each of 
these arrays the two probes used to detect each ORE 
were averaged and treated as single data points with 
p-values determined using a students t-test with a 
one-tailed distribution against the null hypothesis in the 
MATLAB software. Thresholds were set at > 1.5 -fold- 
change for all conditions with the exception of the 
[-Cu/+Cu] and [GE/G] arrays. For those conditions a 
significantly larger proportion of genes were induced, so 
we set the thresholds at > [2a + median] log2 fold- 
change to ensure a similar sized cohort of analyzed 
genes. Genes passing the induction threshold also had to 
pass a p-value threshold of < 0.10. 

Additional files 



Additional file 4: Figure S2. Global Pho7-TAP Enrichment During Pi- 
Starvation. Cells containing the tagged variant of Pho7 (Pho7-TAP) were 
grown in either high-Pi (blue) or starvation (red) media for 120 minutes 
prior to cross-linking and ChlP-Seq processing. Shown are the 
chromosomal enrichment profiles for Pho7 for the S. pombe genome. 
Reads were normalized to total counts for each chromosome. 

Additional file 5: Table S3. Peak List for Pho7-TAP ChlP-Seq Data Set. 
Pho7-TAP peaks detected in both high-Pi (Hi) or no-Pi (No) conditions 
with maximum height > 2x genome average are given. Peak detection 
was performed against a mock sample lacking a TAP epitope as 
described in materials and methods. 

Additional file 6: Figure S3. Pho7-TAP Protein Levels Remain Constant 
During Pi Starvation in cskl^ and csklA Backgrounds, cskl^ or csklA cells 
containing the tagged variant of Pho7 (Pho7-TAP) were grown in high-Pi 
(+Pi) media followed by Pi starvation for 60 or 120 minutes (-Pi). Western 
blot analysis reveals similar levels of the Pho7-TAP protein in all 
conditions as detected by rabbit IgG. The strain lacking the TAP tag is 
shown (Pho7) as a negative control. The Western analysis was completed 
in duplicate, shown is a representative blot. 

Additional file 7: Figure S4. Global Pho7-TAP Enrichment in the 
Absence of Cskl. cskl^ (blue) or csklA (cyan) cells containing the tagged 
variant of Pho7 (Pho7-TAP) were grown in high-Pi media for 120 minutes 
prior to cross-linking and ChlP-Seq processing. Shown are the 
chromosomal enrichment profiles for Pho7 for the S. pombe genome. 
Reads were normalized to total counts for each chromosome. 

Additional file 8: Figure S5. ChlP-Seq Binding Profiles for Pho7-TAP at 
Non-PHO Promoters. Shown are ChlP-Seq profiles for the genes identified 
in Figure 6. As previously described, wild-type cells containing Pho7-TAP 
were grown in either high-Pi (blue) or no-Pi (red) conditions and ChlP- 
Seq libraries were prepared from purified DNA. For comparison, the ChlP- 
Seq signal from mock (black) cells grown in no-Pi and cskl A (cyan) cells 
incubated in high-Pi is included. The gene product of interest is plotted 
based on transcript direction with the plus (+) strand above and the 
minus (-) strand below. Reads were normalized as described in the text 
and the location within the genome is plotted on the x-axis. 

Additional file 9: Table S4. Gene Ontology Enrichment for Pho7 
Dependent Genes Identified via Microarray and ChlP-Seq Analysis. Genes 
possessing a p^oZ^-dependency independent of Pi availability were 
cross-referenced with the Pho7-TAP peak list to identify p/7o7^-regulated 
genes with promoters enriched by Pho7-TAP. Thresholds for induction 
and peak enrichment are described in materials and methods. The gene 
set was processed through the AmiGO toolkit [34,52-55] and GO terms 
containing at least 3 genes with p-values < 0.01 are displayed. 

Additional file 10: Table S5. Microarray Results for p/7o7^-Dependent 
Genes Regulated by Phosphate, Iron, Copper, Osmotic, and/or Carbon 
Switching Stress. phoT^ genes passing the thresholds (described in 
materials and methods) for each stress condition assayed were tabulated 
with the corresponding fold-change (Log2) and p-values (based on at 
least two independent biological replicates). Genes regulated in multiple 
stress conditions by pho7^ are indicated in bold. All comparisons 
between stressed and non-stressed conditions (e.g., [-Pi/+Pi]) were done 
in a pho7+ background. Genes passing the thresholds for each array 
condition have fold-change (Log2) values indicated in bold. +Pi: 10 mM 
H2KP04, -Pi: 0 mM H2KP04, +Fe: 100 uM Fe(lll)Cl3, -Fe: 250 uM DIP, +Cu: 
100 uM Cu(ll)S04, -Cu: 100 uM BCS, 1.2M: 1.2M NaCI, 0.1 M: 0.1 M NaCI, G: 
2% glucose, GE: 2% glycerol, 1% ethanol. 

Additional file 11: Table S6. All primers utilized in this manuscript are 
listed with their primer ID, nucleotide sequence, and experimental 
purpose. 
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Additional file 1: Figure SI. Temporal Dynamics of the Phosphate 
Starvation Response in 5. pombe. (A) Line plot depicting the induction 
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used to delineate the response time are described in the text. Induction 
was normalized to the initial sample pre-starvation (t=0). (B) Average 
expression values for genes in the fast (■) and slow (O) response are 
shown ± SD. 

Additional file 2: Table SI. Microarray Results for Genes Identified in 
the Pi Starvation Time Course. Genes induced following two hours (fast 
response, bold) or four hours (slow response) of Pi depletion were 
tabulated and their normalized expression levels (to t=0 min) are shown. 
Induction thresholds for the fast response were set at > 2o + median 
log2 fold change for each time point (1.00 log2 fold change at 120 
minutes, 1.24 log2 fold change at 240 minutes). Genes above threshold at 
both 120 minutes and 240 minutes post-starvation were classified as the 
rapid response. Genes above only the 240-minute threshold were 
classified as the slow response. 

Additional file 3: Table S2. Microarray Results for Genes Regulated by 
Pi Starvation, pho7^, and/or cskl^. Genes induced following Pi starvation, 
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cells in non-stressed (high-Pi) conditions were tabulated and analyzed for 
orthologs in the 5. cerevisiae PHO pathway. Induction thresholds were set 
at > 1.8 fold-change with a p-value < 0.10. Values shown are the average 
of two independent replicates. 
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